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Abstract. At first, pressure formulas for the electrons under the external potential 
produced by fixed nuclei are derived both in the surface integral and volume integral 
forms concerning an arbitrary volume chosen in the system; the surface integral form 
is described by a pressure tensor consisting of a sum of the kinetic and exchange- 
correlation parts in the density-functional theory, and the volume integral form 
represents the virial theorem with subtraction of the nuclear virial. Secondly on the 
basis of these formulas, the thermodynamical pressure of liquid metals and plasmas 
is represented in the forms of the surface integral and the volume integral including 
the nuclear contribution. From these results, we obtain a virial pressure formula for 
liquid metals, which is more accurate and simpler than the standard representation. 
From the view point of our formulation, some comments are made on pressure formulas 
derived previously and on a definition of pressure widely used. 



PACS numbers: 05.70.Ce,64.30.-ht,52.25.Kn,31.20.Sy 
1. Introduction 

In liquid metals and high density plasmas, their pressure is attributed in main to the 
behaviour of electrons involved in the system and the nuclear contribution is not so 
important. Therefore, almost pressure calculations take account of only the electron 
behaviour in the system; furthermore, the pressure is calculated as that of the electrons 
around one fixed nucleus neutralised in the volume (the Wigner-Seitz sphere). This 
model is called in several ways such as the ion-sphere model, the spherical shell model, 
the Wigner-Seitz model, or the confined atom model. In the case of plasmas, the 
Thomas- Fermi approximation including its extention is applied for this purpose p[]-[p!5[1, 
and in treating metals, the pressure is calculated by solving the wave equation based on 
the density-functional (DF) theory As a consequence, the nuclear contribution 

to the pressure is not taken into account in almost cases. However, even in a high 
density plasma the nuclear contribution becomes important when its temperature is 
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increased. In the present work, we attempt to set up formulas to determine the pressure 
of hquid metals (plasmas) at arbitrary temperature and density by treating them as a 
nucleus-electron mixture using the virial and Hellmann-Feynman theorem. 

There are many investigations |T^-[Q to set up the equation of states (EOS) from 



the nucleus-electron model on the basis of the virial theorem. However, we find confusion 
and improper treatment in some works at the present stage. A typical and important 
example is the virial theorem for the electron pressure Pg in the presence of fixed nuclei 



at the coordinate {Ra} confined in a volume |Tg, which is described as 



SP^n = 2Ts[n]n + E,,{n) - [ n(r)rV^^rfr . (1) 

Here, Ts[n]Q is the kinetic energy defined in the DF theory for the system with a density 
distribution ra(r); Ecs{^) and -ExcH denote the electrostatic and exchange-correlation 
energies, respectively. This equation should be written correctly, as will be discussed 
below, in the form 



SPeil = 2Ts[n]n + E^si^) - / n(r)rV^^rfr 



n 



-5^R„-F,+ i r-P°/-dS, (2) 

which reveals that there are no terms representing the nuclear virial ■ Fq, and a 

surface integral of the exchange-correlation pressure tensor P^^f in the virial theorem 
given by ([I|). The exact equation (Q) involving these terms can be rewritten as 

3P,n = 2 [Ts [n] n + AT^c] + [^es (^^) + Af/,e] -J^RaFa (3) 

a 

= 2Tn + Un-J2Ra-F^, (4) 

a 

where and Uq represent the true kinetic and potential energies, respectively, 
determined by treating the electrons as a many-electron system instead of the DF 
method. Here, the correction terms, ATxc and AUxc, are related to the trace of the 
exchange- correlation pressure tensor Pxc as follows 

2ATxc + AC/xc= [ trPxcrfr. (5) 



n 



On the other hand, the virial theorem ([^) due to the lack of the surface integral term 
leads to the relation |^| 



SE 

2ATxc + At/xc = - / n(r)rV^rfr , (6) 



n 



which becomes zero for a uniform electron density inappropriately. 

The nuclear virial term in (^ causes some confusion in the calculation of EOS, since 
it has an inverse sign as to be considered as the nuclear pressure of the system, or it 



brings a difference from a naive virial theorem. Following Slater [26|, many investigators 
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including works of Janak [|T^], Averill and Painter |]25| and Ziesche, Grafenstein and 



Nielsen I^S] have adopted a definition of the pressure in terms of the nuclear virial 

3P,n = ^ R„-F„ = 2Tn + % . (7) 

a 

The virial theorem (p is a result of this definition of the pressure. In the present work, 
we have proved that the electron pressure in the presence of nuclei fixed at {Rq,} should 
be given by (0) instead of (|^. 

On the other hand, another form of pressure definition is proposed by some authors 
2l| , ^ ^] with use of the Maxwell pressure tensor. However, this definition provides the 



real pressure only when some restriction is imposed, as will be discussed in the present 
work. In the next section, we pick up some important equations to discuss the virial 
theorem for a nucleus-electron mixture. Then, these equations are used to derive the 
electron pressure in the presence of nuclei in Sec. |^, and the thermodynamical pressure 
as a whole system in Sec. Some comments on other pressure formulas are made in 
Sec. § and Sec. ^ The last section is devoted to summary and discussion. 

2. Exact relations to derive pressure formulas 

For the purpose of preparing to derive pressure formulas, we describe in this section 
some necessary and important equations which can be proved without approximation. 

(1) The DP theory: 

Nuclei with the atomic number Z fixed at positions {Rq,} causes an external potential 
f/extlr) for the electrons 

TV 



a=l ' ' ^ ' 

Here, Een is the electron-nucleus part of the total electrostatic energy i^esf^] for this 
system with an electron density distribution ?T,(r) 



2 



Ejn] = — / , ' ^ / drdr' - > / ^—^dr' + - > r (9) 

a=l aitp 



= Eee + Een + -^nn • (10) 

The DP theory provides the effective external potential t/efr(i") in terms of the intrinsic 
interaction part of the free-energy J^i^t = J-'xc + E^e 



Ucs{r) = f/cxt(r) + - - /iint = + - /iint • (11) 

on[rj on{rj on{rj 

This effective external potential provides a way to calculate an exact electron density 
distribution ?7,(r) by solving a wave equation for the one-electron wave function (pi{r) 

n- 



ipi{r) = tiifiir) , (12) 
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in the form 

n{r) = J2fi^^)\^^ir)\\ (13) 

i 

with use of the Fermi distribution /(ej) for arbitrary temperature. 

(2) Electrostatic energy ^'esN- 

Since the interactions among electrons and nuclei are all Coulombic, the electrostatic 
energy -EesN, (iD, is shown to satisfy the following equation 

£^es = - y n(r)rv|^rfr - ^ V^^es (14) 

= - j n(r)rVf/efr(r)rfr + j n(r)rV^rfr - ^ R^V^^es . (15) 

In the above, (|1^) is derived by noting the relation 

rVv{\r - r'l) + r'Vv{\r - r'|) = |r - r>'(|r - r'|) , (16) 

with v'{r) = df(r)/dr, and ( |15|) is led by the use of ([Til). Equations, ([T^) and (|T^), are 
fundamental to derive the virial theorem for a nucleus-electron mixture. 

(3) The kinetic energy in the DP theory: 

In the DP theory, the kinetic energy Ts[n]Q of the electrons confined in an arbitrary 
finite volume fl of the system is defined by the following equation (|l^ in terms of the 
wave functions for (0), and it can be proved that this kinetic energy is associated with 
the external effective potential Ues given by (0) and the density distribution n(r) under 
this external potential in the form 

2T,Mn = 2^/(e.) [ y,*(r)^VV.(r)rfr (17) 



n{r)rVU^s{r)dr + (b A ■ dS . (18) 
n Js 

Here, the tensor A involved in the surface integral for the surface S of the volume fl is 

defined by the following equation (0), which can be represented in other ways such as 

(pop and (21) in terms of the kinetic tensor defined by ( pT]) as follows 

V ■ A ^ ^/(e,)Z^V ■ {^:V[rVv>^] - (rVv9.)(V¥.*)} (19) 
- z^/(e.)^V-MV[rV^,,]-(rV^.)(Vv.*) + c.c.} (20) 

i 

= V ■ [r ■ PT] - I^V^n(r) = V • [r ■ P^^ , (21) 

with 

PZ^ ^ E /(^O^k^V^V.y., + cp.V.VM - (V^¥.n(V.^.) - iV,cp,){V.^:)] . (22) 
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It should be noticed that these relations can be applied to a subspace fl chosen arbitrar- 
ily in the whole space of the system. The proof of the above relations is seen in Appendix. 

(4) Pressure for a noninteracting electron gas in the external potential [/efr(i"): 

The pressure tensor of a noninteracting electron gas in the external potential Ucs{r) 
follows the next equation 

V-P° = -n(r)Vf/efr(r), (23) 

which is a result of the local momentum balance. Here, note that n{r) is the density 
distribution of a noninteracting electron gas in the external potential Uesi^)- 

(5) Mathematical formula: 

For arbitrary tensor P , the following surface integral can be written in a volume integral 
form 



V ■ (r ■ P)dr = f rPdS= / tr Pdr + / r ■ V • Pdr . (24) 
n J s Jn Jn 

3. Electron pressure in the presence of fixed nuclei 

When the nuclei fixed at positions {R^} generate an inhomogeneous electron gas with 
the density distribution n{r), the electron pressure tensor Pe of this system is defined 
by the following equation [|TI], ^ 

N 



V ■ Po = -n(r)V 



-n(r)v|^ . (25) 



-1 ^ 

a=l ' III 

On the basis of this definition for Pg, we can prove the virial theorem for a subspace VL 
arbitrarily chosen in the system as follows 



s 



r-Pe-rfS= <hv{P^' + 'P^^)-dS (26) 



2Ts[n]n + E^,{n)+ [ tr P.^dr - ^ R« ■ F„ 
2Tn + Un - ^Ra-Fa = - j> r • asader ■ dS 



(27) 
(28) 



Q!ef2 

Here, P^^f is the exchange-correlation pressure tensor defined in the DF theory as 

which gives the explicit expression for the electrons in the jellium model at zero 
temperature in the local-density approximation (LDA) in the form 



ide^ 
dn 

with the exchange-correlation energy exc(^) per particle. Also, the force for a-nucleus 



^^^^ = ^^^5.., (30) 



is introduced on the basis of the Hellmann-Feynman theorem: = — Vq-Ees- In (p8|) . 
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the true kinetic and potential energies, Tq = Tg + ATxc and Un = E^si^) + Af/xc, 
which should be provided by treating the many-electron problem, are obtained from the 
correction terms, ATxc and AU^c', these correction terms are determined to satisfy the 
next relations 



2Arxc + AC/xc = / trPxcrfr (31) 
Jn 

ATxc + A[/xc =^xc. (32) 



In the LDA for the electrons in the jellium model, these correction terms are explicitly 
written as 



ATxc = / trPxcrfr-Exc~ / r2(r){3^xc(r) - 4exc(r)}dr , (33) 
Jn Jn 

Af/xc = 2Exc- / trPxct^r^ / n(r){5exc(r) - 3/ixc(r)}rfr (34) 
Jn Jn 

with /ixc = d(nexc)/dn. Furthermore, it is important to note that the virial theorem 
(P5|) can be described in terms of the true kinetic stress tensor ciBader, which is defined 
by Bader for the interacting electrons outside the framework of the DF theory; the 
proof of (^) is given in [^-[^. Thus, the virial theorem with several expressions. 



(^)- (|28|) , provides the electron pressure either in the form of surface integral (|26|) or 
volume integral (^), if the volume is chosen so that the pressure is constant on the 
surface S of Q. At this point, it should be mentioned that the electron pressure in the 
fixed nuclei is represented by the volume integral form (^Tj) correctly if the virial of the 
forces which hold the nuclei is subtracted, while there is no need of this subtraction for 
the surface integral expression (pB]). 

In the following, we give a brief outline of the proof of (p6D -(p8D. At first, we can 
prove that the pressure tensor Pg, (^), of a noninteracting electron gas is identical with 
the kinetic pressure Pk^, (|2T|), in the DF theory 

po = Pf (35) 



with use of (Pl, (gD, (|T|) and (|2T]) as follows 

r-P°-dS= [ trP°c/r+ [ r ■ V ■ P°dr (36) 
Jn Jn 

= 2Ts[n]n - J n(r)rVf/cfr(r)rfr (37) 



r ■ P°^ ■ dS . (3^ 



,DF 

L ■ r 

's 

Next, the electron pressure tensor Pg is shown to be a sum of its kinetic pressure tensor 
P^^ and the exchange-correlation pressure tensor Pxc^ 

Pe = P° + = + PxT , (39) 
from (|23|) and (^) as below 

V ■ (Pc - P°) = n(r) V[f/cfr(r) - 6Ej6n{r)] = n{r)V[6J^,J6n{r)] = V • P°/ . (40) 
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In final, with use of (0), (0), (0) and (|3TD, an identity (^) is sliown to arrive at the 
proof of p6|)-(|28|) as is described in the following 

/ r ■ Pe ■ (iS = / r ■ P° ■ rfS + / r ■ [Pe - P°] ■ = / r ■ (P°^ + P°^) ■ dS (41) 
Js Js Js Js 

= 2TMn - [ n(r)rVt/eff(r)rfr + (f r ■ P^! ■ dS (42) 
Jn Js 

= 2Ts[n]n + EU^)- / n(r)rV-f^rfr 

in ^^(r) 

= 2T,Mn + ^c.(^^)+ / trPxcrfr-5^R„-F„ 



(43) 
(44) 



2[T,[n]n + AT.J + [Eesi^) + AU^,] - ^R„-F„ . (45) 



In the derivation of (|J), Eq. ( poD and the tensor relation (^) are used in the form 

- /"r2(r)rv|^rfr+ / r-P°/-ciS = - f r-V-P°/dr+ / r-P°/-dS = /" trP°/dr,(46) 
J dn{r) Js Jn Js Jn 

which is written in the LDA for the electrons in the jellium model as 

— f n(r)rV-^^^dr = — /" n(r)rVuxc'^r = 3 [n^—e^cdr— (f n'^—e^c^ ■ dS . (47) 
J on{r) J J dn J dn 

4. Maxwell stress tensor and electron pressure tensor 

In the previous section, we have shown that the electron pressure tensor is given by 
Pe = Pk^ + the DF theory. On the other hand, Nielsen and Martin |21| and 



some other [^, |2^ have defined the electron pressure tensor p^'^ using the Maxwell 



pressure tensor as to be 

pNM ^ pDF ^ pDF ^ pM _ ^4g) 

In this section, we inquire into the relationship between these two pressure tensors. 

When the nuclei is fixed at the positions {Rq} with a charge density distribution 
p+(r) = J2a^^^i^ ~ ^a)^ this nuclear distribution produces an electron charge 
distribution P-(r) = —en{r), an electrostatic potential </)(r) of this system is determined 
by the Poisson equation 

V20(r) = -47r[p+(r)+p_(r)]. (49) 

In terms of the electric field E(r) = — V0(r) determined by this Poisson equation, we 
can define the Maxwell stress tensor T as 

T = ( V0 ■ V0 - ^|V0pl ) = -L ( E ■ E - l\E\H] . (50) 



47rV" " 2'"' / 47rV 2 
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Here, 0(r) is represented as a solution of (^) in the form 



M£l±iM£a,,' = _ [ ^ £W (51) 

r — r' r — R„ / r — r' — e 



a=l 



with 



" " ' SE. 



y r-r' ^ r-R„ 

' ' a=l ' 



6n(r) 



(52) 



The Maxwell stress tensor is shown to satisfy the following equation 



V ■ T = -[p+(r) + p4r)]V0(r) = i[p+(r) + p_(r)]V<l>(r) . (53) 



With use of ( |53D and (0), we obtain the relations between the Maxwell stress 

tensor T and the electrostatic energy Ees, (|^), of the system (excluding nuclear self- 
interactions) 

/ r • V ■ Tdr = - / n{r)rV^dr - V R^V^^^es (54) 
Jn Jn on[r) 

= E^{n) = - / tr Tdr + (I r-T-dS. (55) 
Jn Js 

Thus, we obtain the following relation which is important to know how the Maxwell 
stress tensor T plays a role in the virial theorem 

- / trTrfr = — / |Epdr = - / [p+(r)+p_(r)] r ■ V0(r)rfrc/r' - (p r-T -dS (56) 

= E^,{n) - ^ V -T ■ dS. (57) 

At this stage, we introduce the Maxwell pressure tensor = — T; then, we obtain 
the following relations from (|53|), (|5^ and ( pSj ) 

- V ■ (PT + P^! + P^) = 5^ FJir - R„) , (58) 

a 

- / r ■ V ■ (Pf + + P'')dr = J^R^ ■ F„ . (59) 

aen 



Therefore, we get an expression with use of ([57| ) 



/ tr(Pf + P°/ + P^)rfr = 2Tn + Un+ ir-P^-dS, (60) 
Jn Js 

for an arbitrary subspace Q chosen in the system. It should be mentioned that in the 
standard treatment the surface integral in (E7\) is neglected, and as a result, the surface 



integral of the Maxwell pressure tensor does not appear in (|60D; Eq. (^) without the 
surface integral term implies that P^^ defined by (H3) seems to be interpreted as the 



electron pressure tensor giving the virial theorem exactly p2| , p3| . This is not the case, 
since the surface integral exists in (60) in general. At the same time, it should be kept 
in mind that we have the virial theorem for P^'^ 

f r . (Pf + P°/ + P^) ■ rfS = 2T^ + [/f, - 5^ R, ■ F„ . (61) 
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only if a volume Q is chosen so as to fulfill the condition r • • dS = on the surface 
of this volume; in this case, ( |6TD reduces to (p^). 

In conclusion, as a consequence of (|60|) and (|6TD, the tensor P^^ defined by (^) 
plays a role of the electron pressure tensor in the presence of the fixed nuclei only for a 
subspace fl, on which surface the condition r ■ P^ ■ (iS = is satisfied. 



5. Comparison with other electron-pressure formulas 



In Sec. y, we derived the virial theorem (^61)- (Pq), in the surface integral and volume 
integral forms, for the electrons in the presence of fixed nuclei, which provides the 
electron pressure in condensed matter at arbitrary temperature and density. On the 
other hand, there have been proposed many methods to calculate the electron-pressure 
in the fixed nuclei. In the present section, we will make comparison of our result with 
some of other formulations, and give some comments. 

If the electron pressure is assumed to be constant on the surface of the Wigner-Seitz 
sphere, Eq. (^) with (^) provides the electron pressure in the presence of the fixed 
nuclei in a surface integral form 



r ■ P, ■ rfS 



3P,V 



{v7*V[rV<^i] - (rVv?i)(V¥P*) + c.c.]dS 



+ 



r P 



DF 



dS. 



(62) 



When we use the LDA expression for P^^^^, — n^dexc/dn^^j, in the above, (^) reduces 
to Liberman's pressure formula. However, the nuclear virial term does not appear in 
(|5^) on comparison with the Liberman expression (7). 

On the other hand, if we choose the volume VL to be the Wigner-Seitz sphere with a 
radius Tq, the virial theorem in the surface integral form, (|26|) , with (|22|) can be rewritten 
as 



P. 



s 



cff 



2mrQ 



( ^^^^ 



fi_ 
ro 



3? 



dVt 
An 



+ P°/(ro) (63) 



with dQ = sin 6d6(f) and the definition of an operator £'^Y^{6, (p 

Equation (^) becomes the pressure formula given by Nieminen- Hodge |]18[ and Pettifor 
|19| . Furthermore, when the volume is chosen to be the Wigner-Seitz cell with 
the muffin-tin approximation, and the exchange-correlation effect is described by a 
generalized gradient approximation, Eq. (|6^) leads to the pressure formula given by 
Legrand and Perrot by taking account of the discontinuity in the electron density 
distribution. 
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Nielsen and Martin with other investigators |2^, have defined the pressure 
in terms of the Maxwell pressure tensor in the form 

pNM^pDF^pDF^pM_ ^g4) 

As was discussed in details in Sec. ^, this tensor P^'^ cannot be interpreted as the 
electron pressure except the case when we treat the pressure in a restricted subspace fl 
chosen to satisfy the condition ^^r ■ ■ dS = on its surface (e.g. Wigner-Seitz cell). 
When this condition is fulfilled, we obtain virial expressions such as 



trP^^dr 



r P 



NM 



dS 



2Tn + Un , 
2Tn + Un - 



(65) 
(66) 



In addition, Ziesche, Grafenstein and Nielsen have defined the pressure in terms of 
R^-Fa considering a solid as an isolated system [that is, r ■ P^m ■ rfS = 0] by 

3P,n = Y,^a-Fa = 2Tn + Un, (67) 



as was defined by Slater |^. However, the electron pressure in the fixed nuclei should 
be defined by (pSf), where the virial of the forces that hold the nuclei fixed is subtracted 
here to be added later in the thermodynamic pressure formula ( pB]) as will be discussed 
in the next section. 



Janak [0 derived the pressure formula starting from the virial theorem of the 
volume integral form which is written in our version as 



3PV = 2Ts[n] + ^es + / trP^Jr . 

Jv 



{61 



In this equation, he replaced E^s in a volume V containing nuclei by an approximate 
expression (|7lD 

(69) 



Ecs 

— / n(r)rV?7efT(r)(ir - 

V 

N - n(r)rVf/eff(r)dr 



n[r)r\/———dr 
dn{r) 



Ra Va-Ee 



n(r)rv|%dr - 



(70) 
(71) 



1 z-^ 

as was written in his equation (18) for the unit cell fio- Thus, the Ewald term —-^C-^ 
can be taken as the nuclear virial in this approximation. With this replacement, Eq. (^) 
is rewritten as the next equation ([72D, which leads to a simplified expression of his 
pressure formula [his (24)] as (|73D 



3PV 



r-Pe-(iS + ^Ra-Fc 



3A^ 



2 n , C^Cxc 

-ep + n— — 
5 dn 



1 72 

4 a 



(72) 
(73) 
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with ep being the Fermi energy of a noninteracting electron gas, if the muffin-tin 
approximation is not introduced. In the derivation of (|73D, the electron density n is 
assumed to be constant at the surface of the unit cell Qq and we used the following 
relations 

Ts[n]n, - TsKU = \ I n{v)TVU,^{v)dv , (74) 

and 

Ts[n,]no = \j nb(r)rVf/efr(r)rfr, (75) 

for the bound-electron density nb(r). Equation (^) indicates an essential structure of 
his pressure formula, although simplified. For the same reason, his formula [his (25)] for 
the total energy becomes 

N f \ 7'^ 

E = -- n(r)rVf/eff(r)dr + E[nl] - {AT^Mr)] - AT^cK]} - N-C^ , (76) 
^ Jfio 4 a 

where ATxc[^(r)] in the LDA is given by (p^). 



6. Pressure in liquid metals and plasmas 

In Sec. 1^, we have derived pressure formulas for the electrons in the fixed nuclei. 
However, the thermodynamical pressure of liquid metals and plasmas must be 
determined as a nucleus-electron mixture as a whole; the electron pressure Pe is only a 
part of the whole pressure P, although Pg is nearly equal to the total pressure in metals. 

In general, the thermodynamic pressure is defined by P = —dF/dQ\N,T using 
the free-energy of the system: F = —ksT In Zn with the partition function Zn = 
Tt exp{—f3H) of a nucleus-electron mixture. Here, its Hamiltonian H = Hee + Hen + H^^i 
consists of the electron and nucleus Hamiltonians and -ffnn, in conjunction with 
the electron- nucleus interaction Hamiltonian H^^. When the nuclei can be treated as 
classical particles [Pq,,Rq,] = 0, the partition function is written in the form (see 
for details) 



jdR^J dP^ exp[-pH^s] 



(77) 



with 



a 

<I>({Rj;n) = Jo + ^es[ni,n] +^,Jni,n] = j;[{Rj;n]. (79) 

Here, ni(r) = J2a ~ -^o) denotes the microscopic nuclear density distribution. Also, 
it is important to realise that $({Rq,}; n) plays a role of a many-body interaction among 
the nuclei, and jFe[{RQ,};n] = ^({Rq,};^) represents at the same time the free energy 
of the electrons under the external potential caused by the fixed nuclei at the positions 
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{R„} in the DF theory with reference to jFg, the free-energy of noninteracting electron 
gas. 

A volume change of the free energy can be calculated by introducing a new variable 
A defined by f2(A) = X^Q in the following way 

Aap 

3~dX 



dn~ ^ ' dx dn{x) 



A=l 



(80) 



A=l 



Here, F{X) = —k^T In Zn{X) is defined by 

ZN{n) = j^JdR''J rfP^exp[-/5i/,ff(R,P)] 



n{x) 



N 



dR^ I dP^ exp[-pH^siv^,P)] = ^iv(A) 



^1) 

J^lfl^N I "'^^ / ^^^^L A---onv/--7- /J — -iVV'V v32) 

with R = R/r7(A) and ri^{X) = Q{X). Then, the A-derivative of F is easily performed as 
follows 



dFjX) 
~dX~ 



Q 'dX 



n{x) 



N 



dR^ I dP'^ exp[-l3H^s{7]R, P)] 



with 



eflf 



dX 



Zn{X) Jq 
9$(r/R,fi(A)) 



cS 



dX 



A=l 



dX 



A=l 



J]R„V„<l' + 3^]— . 



(83) 



^4) 



Finally as a consequence of the above equation, the thermodynamical pressure of liquid 
(solid) metals and plasmas is established as a mixture of nuclei and electrons in the 
form 



3PQ = -3Q 



dF 



- X 



dF 
dX 



A=l 



3NkBT + ^ I c/R^exp[-/5<l>] i -3fi|| + ^R„-F, 
3NkBT+(-3n^) + ( ^R„-F, 



with 



1 



rfR^exp[-/?$({Rj)], 



(85) 
(86) 

(87) 



and A being the de Broglie thermal wavelength of the nuclei. Here, —VLd^^/dVL = P^Q is 
the electron-pressure in the fixed nuclei, which can be determined by ( p6D or (^) within 
the framework of the DF theory, since $ = jFe[{RQ,}] is the free energy of electrons in 
the presence of nuclei fixed at {Rq}. At this stage, an important remark should be made 

11 that the virial theorem for 



(89) 



concerning (|86D . It has long been recognised |12, ^ 
a liquid (or solid) metal is written in the form 

3-Pc^ = Ra'Fo = 2Tq + Uq . 
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as was done by Slater [^. However, Eq. ( pB]) indicates that the electron pressure under 
the external potential caused by fixed nuclei should be given by 



3Pen = -31^^ = 2Tn + Un-J2^--^-^ 



(90) 



aefl 



as was shown in (pH]), and that the thermodynamical pressure P of a liquid (solid) metal 
should be defined by (^) as a nucleus-electron mixture, where the nuclear virial term 
in ( ppp disappears by cancellation. 

Thus, the EOS of the nucleus-electron mixture is determined by two types of the 
formulas using a surface integral and a volume integral in the final form 



3Pn = SNkBT + jr- (P^ + PD/^ ■ dS + R„-F„ 
= 3iVfcBT + (2T,[n] +E,, + j trP°/dr' 



(91) 
(92) 



with combined use of (|86D , (|26|) and (pTl) . It is worth pointing out that the volume 
integral form ( p2D of the virial theorem contains the nuclear virial contribution, although 



it does not appear explicitly comparing with the surface integral form (0). In addition, 
it should be kept in mind here that the term P^^f (ni, n) depends on the ion- configuration 
ni(r), and therefore it is not that of the jellium model in a strict sense. In the 
jellium model with use of the LDA, the exchange-correlation pressure is given by 
P^Iav = n^de^d dn5^y. In this approximation, we obtain the following expressions 



trP^/rfr 



n 



dn 



pDF , p 



df\ 

xc / 



■dS^3n 



2 I 2 '^^xc 

n-ep + n — — 
5 dn 



(93) 
(94) 



where the electron density is assumed to be constant n^ on the surface of the volume Q. 

On the other hand, in liquid metals and plasmas where their effective interactions 
among ions with uniform density nj^ can be described by a pair interaction vff (r), the 
average of the nuclear virial is given by 



NnQ I gu{r)rVvjf{r)dr 



(95) 



Therefore, the EOS of liquid metals and plasmas is described in this approximation by 
a formula 



If 1 

+ g^o / g^ceir)vec{r)dr - -{nlf 



"f" ^xc 



^ii(r)rVfif (r)rfr (96) 
(7n(r)rVt;if (r)rfr. (97) 
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Here, g'^^.ir) is a free electron part of the electron-electron radial distribution function 
jy] defined by 



l + n^o^Q[gl{r)-l] 



— dio 



hp 



{9t 



in terms of the dynamical local-field correction G{Q,u!) in the electron-ion mixture. In 
the jellium approximation, we obtain expressions, G{Q,u!) ~ G^^^^{Q,u!) and 

'To 



n 



jell 



3 I o '^^xc 



n 



(99) 



to be used in (^) leading to (pTf). At this point, it should be noted that when we 
approximate $({Rq,}) at earlier stage to be 



we obtain another EOS 



idUp 
' drin 



6^ 



rV — 3n, 



I d 



dni 



v!^{r)dv. 



(100) 



(101) 



which is used as a standard virial formula for a liquid metal |41]. The difference between 
(p7|) and (|101|) comes from the electron-pressure term in (|87|). This term with an 
approximate expression to $ ( [LOOj) leads to (|101|) , while this term with use of the exact 
expression (p6D yields a more accurate and simpler pressure-formula 



7. Summary and discussion 

As a first step, we have derived pressure formulas, ( p6D and (p7|), for the electrons under 
the external potential produced by fixed nuclei on the basis of the DF theory; one is 
a surface integral form and the other, a volume integral form showing the exact virial 
theorem (^). The electron pressure tensor is represented by a sum of the kinetic tensor 
and exchange-correlation pressure tensor, Pg = P^^ + Pxc^, in the DF theory. 

As a second step, the thermodynamical pressure P is defined by the volume 
derivative of the free energy for a nucleus-electron mixture (liquid metals or plasmas) 
on the assumption that the nuclei behave as classical particles. This procedure leads 
to the EOS formula given by (^61), which contains the electron-pressure in the presence 
of fixed nuclei and the nuclear virial. In final, we obtain pressure formulas, (^Tj) and 
(^), for liquid metals and plasmas in the surface integral and volume integral forms, 
respectively. In this respect, it is important to recognise that the exchange-correlation 
pressure tensor P^^{ni,n) is that of the electron- ion mixture, which is approximated 
by that of the jellium model in the standard treatment. 

For a liquid metal, when we can approximate a many-body interaction $({Rq,}; n) 
([f9|) among the nuclei to be represented by the sum of a binary interaction wff (r), 
Eq. (|9ll) is reduced to a virial pressure formula (^). Furthermore, when the electrons 
in the mixture can be treated by the jellium model, (that is, the electron is weakly 
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coupled with the ions in the electron-ion mixture), we obtain a virial pressure formula 
(p7|), which is simpler and more accurate than the standard virial formula ( |101| ) for a 
liquid metal . 

Also, our formulation for the pressure have clarified confusion and improper 
expressions found in the previous works, as was addressed in the introduction. The 
electron pressure in the presence of fixed nuclei should be determined by ( pll ) instead 
of dH), as was discussed in Sec. |. Also, a tensor, P^^ = pDF ^ pOF ^ pM ^^^-^ ^^le 
Maxwell pressure tensor P^, can be interpreted as a pressure tensor only when we treat 
a subspace Q of the system fulfilling the condition r ■ P^ ■ rfS = on its surface. 
Therefore, there is no need to introduce a correction term due to the Maxwell pressure 



tensor in the Liberman pressure expression |jT6[, in spite of the suggestion by FolUand 
2^ and Nielsen and Martin . 



There is a fundamental problem whether we can choose arbitrarily any subspace 
in the system for the virial theorem to be applicable. Srebrenik and Bader p9| have 



answered to this question by giving the zero-fiux surface condition to a subspace; a 
subspace has a close surface through which the flux of Vn(r) is everywhere zero, that 
is, 

Vn(r) ■ e(r) = 0, (102) 

for any point r in a close surface with a normal vector e(r). In a local virial theorem 
constructed on this condition, there is no ambiguity in the definition of kinetic energy 
Tg and the kinetic tensor P^^, as will be discussed in Appendix. 
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In the DF theory, the wave function is determined by the Schrodinger equation with a 
self-consistent external potential t/efr(r) 



2m 



V' + f/eff(r) - e. 



0. 



(A.i; 



By multiplying (|A.1| ) by V5*rV, there results 



2m 



2m 

+ [U,s ~ ei]ip*rVipi = 0. 



(A.2) 



Therefore, from the above equation we obtain 



2m 



[2^*VV.-V{(vpn'V(rv..M*)}] 



This equation is rewritten in the form 



with 



A 



2m 
Am 
Am 
Am 



2m 

V-M*V[rV^.]-(rV(/.0(V^:)} 

V ■ {^*V[rVvP,] - (rV^,)(V^:) + c.c.} 

V ■ {V(^:^,) + (rV)V(^*^.) - 2[(rVvP,)(V^*) + (rV^*)(V^.)]} 

V ■ {V(^:^,) + ^*(rV)V^, + ^.(rV)V^* 



(A.3) 
(A.4) 

(A.5) 
(A.6) 
(A.7) 

(A.8) 



-(rV^,)(V^*)-(rV^n(V¥^.)}, 

which leads to (|l|) with (|19D-(|l]). 

On the other hand, (p!8D can be proved in a direct way by following Pauli [43|. Pauli 
have proved that the probability current 



. _ 1 

^ " 2 



m m 



(A.9) 



for a wave function of a single electron moving under the external potential t/efr(r) obeys 
the following continuity equation 



d 

m—i = -V ■ S - {ipip*)VU,s{r) 



(AlO) 



with a tensor defined by 



-re 



Am 



bV^V,<^* + <^*V^V,<^ - V^^V.v* - . (A.ll) 
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This provides the following equation for the static case 
d 

mr ■ — i = — r ■ V ■ S — {ipip*)rVUcs{r) = . 
When we interpret defined by 



as a kinetic tensor in the DF theory, (|A.12| ) is changed into the form 

n{r)\/Ucs{'r)dr 



r ■ V ■ P^^dr 



trP^'^rfr - ^ r ■ P'^' ■ dS 
with use of a tensor relation (^). By noting the relation 



DF 



trS = J2 "^MM 



Am 



2m 2m 



4m 



2m 



2m 2m 



1 




2m 


2 



we can represent the trace of P^^ in the form 

[ trP^^'rfr = 2TAn]n + — i Vn(r) • dS . 
Jn 4m Js 

From ( |A.14| ) and ( [A.19| ), we obtain in final 



+ 



n{r)VU,s{r)dr 



r • P°^ ■ dS 



4m 



Vn(r) ■ dS , 



which is identical with (|1^) with ([2T|). In this respect, it is interesting to 
have arbitrariness to define the kinetic-energy functional in another form 



2Ts[n] 



trP°^rfr 



Jn 



dr 



+ 



4m 



V^?T.(r)(ir . 



involving the last term of (|A.20|) here. 
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